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ABSTRACT 

Using broadband optical imaging and Chandra X-ray data for a sample of 46 cluster central domi- 
nant galaxies (CDGs), we investigate the connection between star formation, the intracluster medium 
(ICM) , and the central active galactic nucleus (AGN) . We report the discovery of a remarkably sharp 
threshold for the onset of star formation that occurs when the central cooling time of the hot atmo- 
sphere falls below '--^ 5 x 10^ yr, or equivalently when the central entropy falls below ~ 30 keV cm^. In 
addition to this criterion, star formation in cooling flows also appears to require that the X-ray and 
galaxy centroids lie within ~ 20 kpc of each other, and that the jet (cavity) power is smaller than the 
X-ray cooling luminosity. These three criteria, together with the high ratio of cooling time to AGN 
outburst (cavity) age across our sample, directly link the presence of star formation and AGN activity 
in CDGs to cooling instabilities in the intracluster plasma. Our results provide compelling evidence 
that AGN feedback into the hot ICM is largely responsible for regulating cooling and star formation 
in the cores of clusters, leading to the significant growth of supermassive black holes in CDGs at late 
times. 

Subject headings: galaxies: clusters: general — galaxies: elliptical and lenticular, cD — X-rays: 
galaxies: clusters — cooling flows 



1. INTRODUCTION 



X-ray observations of the intracluster medium (ICM) 
in a majority of clusters show stongly pe aked central 
emission from gas with short cooling times (Peres et al. 
1998 1. Such c lusters were thought to harbor cooling flows 



( |i:''abian|1994[ |, in which gas cooling out of the ICM at the 
core ot the cluster is replaced by gas at larger radii in a 
slow inward flow. Low-spectral-resolution X-ray observa- 
tions predicted_de£Osition rates of 100— lOOO's M© yr~^ 
(e.g., Peres et al.|1998 ). These rates and the implied cold 
as masses exceeded the inferred star formation rates 



McNamara 19971 and the observed gas masses (Edge 
UUip in the central galaxy by an order of magnitude 



i'his apparent disagreement between the amount of gas 
seen cooling and the amount seen in traditional sinks 
persisted for more than two decades. 

Recently, however, high-spectral resolution observa- 
tions by the Chandra and XMM-Newton X-ray Obser- 
vatories have shown that the cooling rates predic ted by 
ROSAT are too large by a factor of ten (e.g., David 
et al.||2001| [Blanton et al.||2003| [Peterson et al.||!^00^| 



Kaastra et al. 20041. These new telescopes have shown 
that cooling-flow spectra lack the emission signatures ex- 
pected from gas cooling below ~ 1 keV. Heating by the 
central active galactic nucleus (AGN) has emerged as the 
most promising means to prevent signiflcant amounts of 



cooling from occurring (for a review, see McNamara & 
Nulsen|[2007t . 

Some form of heating is also required in hierarchical 
CDM simulations to reproduce the high-luminosity end 
of th e galaxy luminosity function . Recent simulations 



(e.g., |Croton et al.||2006| |De Lucia fc Blaizot|[2007| have 



found that strong feedback at the cores of the largest 
halos is necessary to prevent over-cooling and the forma- 
tion of much more massive galaxies than are observed. 
The local galaxy luminosity function show s an exponen- 
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tial cutoff at the highes t luminosities (e.g., Benson et al. 
Bower et al.l2006 1 that is reproduced well by mod- 
b 



hierarchical formation that include AGN feedback 
However, some net cooling is probably occurring, as 
evidence of active star forma tion is observed at the c ores 
of many cooling flows (e.g. IJohnstone et al.||1987j |Ro 
manishini 1987: Allen 1995; Smith et al. 199?'; HicksT2 
Mushotzky|2005i|Donahue et al.|2007f|O'Dea et al.|2008| 



and it has been known tor some time that indicators of 
star formation in the centrally dominant galaxy (CDC), 
such as spatially extended excess blue emission and opti- 
cal line emission, correlate, alb eit roughly, with the prop- 
erties of the cooling flow (e.g., | He ckmanlll981| iHu et alT] 
l985| [McNamara fc O'Connell|ra 89 Cardiel et al. 199^ 
Additionally, CDGs with strong optical line emission are 
preferentially found near the cores of cooling flows, im- 
plying a physical connection between star formation and 
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the c ooling ICM (Crawford et al. 1999 Edwards et al. 
20071 



However, only recently, with the advent of Chandra, 
has it become possible to derive the X-ray properties on 
the same scales as the optical properties to determine 
more precisely the relationship between cooling and star 
formation. In this paper, we use optical and X-ray data 
of similar spatial resolution for 46 CDGs to examine the 
connection between the ICM properties, such as the cen- 
tral cooling time and the AGN's feedback power, and 
the presence of star formation in the central galaxy. We 
adopt Hq = 70 km s'^ Mpc"\ f^A = 0.7, and Q.m = 0.3 
throughout this paper. 

2. THE SAMPLE 

The sample comprises 39 CDGs for which we have ob- 
tained broadband optical data, plus an additional 7 sys- 
tems for which we have taken optical data from the lit- 
erature. The sample was constructed to test whether 
the presence of active star formation in the CDG is re- 
lated to the properties of the cooling flow. To that end, 
the objects in the sample were chosen from the Chandra 
Data Archive to include both cooling flows and noncool- 
ing flows, with a wide range of central cooling times, from 
icooi ^ 10^ yr to tcooi ^ 10^" yr. Additional constraints 
on the sample were that the available X-ray images have 
sufficient counts to derive density and temperature pro- 
files and that the systems be visible to the optical ob- 
servatory during the alloted observing windows. Table [l] 
lists the general properties of the objects in the sample. 

3. OBSERVATIONS AND DATA REDUCTION 

3.1. Optical Data 

The primary signatures of star formation are due to 
the presence of hot, massive stars in the star-forming 
population. This population is mixed with the generally 
old stellar population of the underlying CDG, as well as 
with cold gas and dust. These massive stars emit most 
of their energy in the ultraviolet and blue parts of the 
spectrum, unlike the less massive but much longer-lived 
stars that account for most of the CDG's emission at 
longer wavelengths. The ultraviolet emission in turn can 
ionize nearby cool gas, res ulting in optical line emission 



such as Ha emission 



Heckman et al. 1989 Voit & 



Donahue|p97l [Hatch et al. ||2007 p . Additionally, the en 
ergetic photons from the massive stars heat surrounding 
dust , which re-emit s the light at far infrared wavelengths 
( Q'Dea et al.|2008" I . In this study, we use the excess blue 
emission (above that expected from the old background 
population) as an indicator of recent star formation. 

We use broadband imaging at short and long wave- 
lengths to search for the excess blue emission in the 
cores of the CDGs. By comparing the short-wavelength 
surface-brightness profile, which is sensitive to young, hot 
stars, with the long-wavelength profile, which traces the 
old background population, we can detect the presence 
of star formation. Elliptical galaxies generally show no 
signs of recent star formation, and have color profiles that 
become bluer with increasing radius (e.g., [Vader et al. 
TW; 'Fra nx et al.][T989] [Peletier et ar||1990| [Goudfrooij 
et al. 1994|. This effect is thought to be generally due 
to a decreasin g stellar me talicity wi th increasing radius 



color profile that becomes increasingly blue towards the 
center is unusual (i.e. the CDG has excess blue emission), 
and is indicative of recent star formation in that region. 
We use this property to identify star formation in our 
sample. 

The excess blue emission, while generally indicative of 
star formation, can be the result of other sources, such 
as the scattered light from an AGN or a low-metalicity 
stellar population. However, because AGNs are point 
sources, their emission should not be spatially extended 
on scales of tens of kpc, as the excess blue emission is 
often observed to be in CDGs. Scattered light can make 
the AGN emission appear more extended, but studies of 
the nearby cooling fiows A1795 and A2597 have ruled out 
significant polarization in the excess blue light (e.g.. Mo 



Namara et al. 1996a 1999 ) . As to metalicity effects, most 
modeling of spectra and broadband colors of CDGs has 
found that systems with large color excesses are better 
modeled by emission from massive O and B stars, rather 



than old, low-metalicity stellar p opulatio ns (e.g., Allen 
|T995l |Cardiel et al.|199ll|Smith et al. 1997l|Cardiel et a 



1998 ; Crawford et al. 111999 1. Additionally, Hubble Space 



Telescope studies of two CDGs with prominent blue ex- 
cesses, A1795 and A2597, have resolved the blue emission 
into knots and demonstr ated that it is mos t likely due 
to recent star form ation (Koekemoer et al. 1999, ,2002] 
Q'Dea et al.|[2004| . 



e.g., [Carollo et al. 1993[ Kobayashi fc Arimoto 



1999 



TamurrSrOhta|2006| |l''erreras et al.|2005p . Therefore, a 



Since our goal is to search for connections between star 
formation and the cooling flow, we used data that are 
sensitive to the presence of active star formation. Such 
star formation is best detected at short wavelengths due 
to the rising spectral energy distribution (SED) of hot 
stars in the ultraviolet and the falling SED of the back- 
ground galaxy population. Therefore, we chose to ob- 
serve in the U band (with a central wavelength of 3582 
A) . This choice minimizes contamination from redshifted 
[OIIJA3227 emission lines expected to be present in many 
of the systems. We also need long- wavelength observa- 
tions to trace the old stellar population; however, at long 
wavelengths, contamination from Ha and H/3 lines can be 
a problem; therefore, we chose to observe in both R and 
/ fllters (with central wavelengths of 6513 A and 8204 A, 
respectively). Colors obtained from both fllters give us 
a means of verifying that any anomalous blue emission 
is due to star formation. If the U — R and U — I colors 
are both anomalously blue, it is likely to be due to star 
formation, and not due to contamination from emission 
lines (at least in the red images). 

Optical data were obtained during flve separate runs. 
Three runs, totaling 11 nights, were done using the 2.4- 
m telescope at the MDM observatory on Kitt Peak, Ari- 
zona between the dates of March 9-12, 2005, September 
26-29, 2005, and May 22-24, 2006. An additional run 
of four nights was performed at the KPNO Mayall 4-m 
telescope between October 5-8, 2005. Lastly, a 2-night 
run was performed at the WIYN 3.5-m telescope over 
January 22-23, 2006. See Table [2] for details of the opti- 
cal observations. Broadband imaging was done at MDM 
using the Echelle CCD with a 9.4 x 9.4 arcmin fleld of 
view and a scale of 0.28 arcsec pixel^^. At the KPNO 
4-m telescope, imaging was done using the Mosaic CCD 
array of eight CCDs with a 36 x 36 arcmin flcld of view 
and a scale of 0.26 arcsec pixel"^. At the WIYN 3.5-m 
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TABLE 1 
Sample Properties. 
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Note. — The X-ray and CDG core positions are from this work; see Sections |4.1| and |4. 2. l| Systems that lack optical CDG positions were 
not observed in this study; optical data for these systems was taken from the literature (sec references given in Table jsj. Uncertain core 
positions (due to the lack of a cooling core or to the presence of confusing substructure) are marked with an asterisk. 
^ Total X-band magnitudes from the 2MASS catalog. 

^ Projected radius from the X-ray core to the CDG's core. Separations of less than 2 aresec are treated as upper limits. 



telescope, imaging was done using the Mini-Mosaic CCD 
array of two CCDs, with a 9.6 x 9.6 arcniin field of view 
and a scale of 0.14 arcsec pixel"^. 

Harris U, R, and / filters were used during all runs. All 
objects were imaged in the U band and in either i? or / 
band (most objects were imaged in both R and / bands). 
Exposure times were typically 2 x 600 s in J7, 2 x 150 s 
in /, and 2 x 225 s in R; however, in some cases (e.g., 
for distant systems), longer total exposure times were 
required to achieve sufficient signal to noise. Multiple 



frames were taken to allow for easy removal of cosmic- 
ray events. The frames were dithered by 30 — 60 arcsec 
between exposures to allow for the removal of artifacts 
such as bad columns and CCD gaps and to improve flat 
fielding. Conditions varied during the runs, with pho- 
tometric conditions for approximately two-thirds of the 
total time. 

Optical reduction and analysis was performed using 
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TABLE 2 
X-RAY AND Optical Observations. 
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IRAF (the Image Reduction and Analysis Facility^), ver- 
sion 12.2.1, and custom procedures written in IDL.^ 

3.1.1. Standard Reductions 

The bias level and its row-to-row variation were mod- 
eled and subtracted in all frames by fitting a high-order 
function to the overscan region of each chip. The re- 
maining bias structure was removed with a bias frame 



1 See 

2 See 



http://iraf.noao ■ edu/ 



http : //www . ittvis . com/idl/ 



constructed by averaging 30-40 individual frames taken 
each night in the evening and morning. Twilight flats 
were then used to remove differences in pixel sensitiv- 
ities across the chip. Flats were made for each night 
by averaging together at least 3 dithered frames taken 
during the evening and morning twilight with exposure 
levels of ~ 1/2 the saturation level of the CCDs. 

The MDM Echelle CCD was found to have signifi- 
cant dark current 0.1 — 0.2 counts s^^) during the 
March and September 2005 runs (the dark current dur- 
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ing the May 2006 run was much lower) , which resulted in 
large sky gradients in the long-exposure J7-band frames. 
Additionally, a small but significant light contamination 
(probably due to light from the red light-emitting diodes 
on the instrument mounting) was present in the R- and 
/-band dark frames. We corrected for the dark current 
and light leak by creating a master dark frame in each 
band by averaging 4-5 individual frames together and 
subtracting a scaled version of this master frame from 
each object frame. Each individual dark frame had a 
dark time typical of the longest exposure time used dur- 
ing observation (e.g., 900 s for U). Unfortunately, the 
dark current was found to vary both in magnitude (by a 
factor of two) and in spatial structure over the course of 
a single night. This variation made it impossible to use 
a simple scaling by dark time. Instead, the best scaling 
was determined using a minimization routine. 

The best dark scaling for each object frame was deter- 
mined by minimizing the residuals between the modes 
of the pixel values in multiple regions of the CCD af- 
ter dark subtraction. This method effectively scales the 
dark frame to produce the flattest sky across the CCD. 
However, for this method to be effective, sky variations 
due to pixel-to-pixel differences in sensitivity must first 
be removed by flat fielding. Therefore, the dark images 
were also fiat fielded before fitting and subtraction. The 
fitting was done in IDL using the MPFIT package.'^ Typ- 
ically 9 object-free sub-regions of the chip in which the 
dark structure was most pronounced were used for fit- 
ting. If the fitted scaling in any given region differed 
by more than 3 sigma from the adopted mean it was re- 
jected. The mean scaling across all the remaining regions 
was then used for dark subtraction. Remaining sky gra- 
dients were typically on the order of 2-3% of sky across 
the CCD, but they are generally much smaller 0.5%) 
in the region of the BCG. 

Fringing due to night sky lines was found to be sig- 
nificant in the /-band images. Therefore, a fringe frame 
for each night was constructed and subtracted from the 
/-band object frames. The fringe frame was made by av- 
eraging together all the /-band object frames of a given 
night, after masking of all objects (using the objmasks 
task in IRAF). The sky level and any large-scale gradi- 
ent were removed from the fringe frame using a median 
filter over large blocks of pixels to filter out small-scale 
variations. The final fringe frame was then scaled and 
subtracted from all the /-band object frames using the 
rmfringe tool in IRAF. 

In the KPNO data, an additional additive feature was 
present in /- and J7-band images due to scattered light off 
of the prime-focus corrector (this feature is not present 
in the /?-band images due to the anti-refiective coating 
applied to the corrector) . A template image of this "pupil 
ghost" was made using the mscpupil task in IRAF. This 
task isolates the pupil contribution to an image by fitting 
a spline function to the background of a master image 
and subtracting it off. The master image was made by 
averaging all flat-corrected images in a given filter over 
the entire run, after the masking of all objects using the 
ohjmask task in IRAF, and scaling by the mode of pixel 
values in object-free regions of the image. The resulting 
pupil template image was then scaled and subtracted off 



of each image interactively. 

Once the additive features were subtracted from each 
frame, a dark-sky flat was made by again averaging all 
frames in a given filter over each night. Generally, 10- 
20 frames were required to produce good dark-sky fiats. 
When an insufficient number of frames existed for a given 
night (due, e.g., to poor weather), a dark-sky fiat from a 
neighboring night was used. Once again, all objects were 
masked before averaging. The sky fiats were then median 
smoothed with a 129 x 129-pixel box filter to reduce the 
noise but preserve large-scale structure. 

Finally, a world coordinate system (WCS) solution was 
derived for each frame using the WCS mapping tools in 
IRAF. The WCS was accurately aligned and scaled to 
each frame using fits to USNO A2 catalog positions of 
typically 100-200 stars. After this adjustment, the WCS 
was found to be accurate to ±0.3 arcsec or better. 

3.1.2. Sky subtraction 

Accurate determination of the sky surface brightness is 
essential for tracing galaxy profiles to large radius, where 
the galaxy's surface brightness is often just a few per- 
cent of the sky. Our study is concerned mainly with the 
profiles in the centers of the galaxies, where their emis- 
sion dominates over the sky. However, we have nonethe- 
less attempted to determine the sky level as accurately 
as possible. To this end, we have adopted the method 
of sky determination used in |McNamara fc O'Connell 
(19921, in which counts in source-free regions of the sky 



are assumed to be dominated by Poisson noise. The dis- 
tribution of counts in each region should then be well 
represented by a gaussian whose width is an estimate of 
the statistical error in the sky level. Therefore, to deter- 
mine the sky level, we extract the counts in a number of 
source-free regions, construct a histogram of sky values, 
and fit a gaussian to the histogram. This method of fit- 
ting the distribution of sky values, in contrast to using 
a simple mode or mean, allows one to use the distribu- 
tion's deviation from a normal distribution to identify 
the presence of contaminating sources. 

We use the mean of the best-fit gaussian as the modal 
sky value for a given region. We used 7-9 sky regions 
with areas of at least 1000 pixels each, distributed over 
the chip in the region of the CDG (but far enough away to 
minimize contamination), and eliminated regions whose 
sky modes were more than 2cr from the adopted mean. 
Comparison with sky values determined by fitting a plane 
to the entire CCD (median filtered over large blocks) 
showed typical differences in sky values between the two 
methods of < 1 — 2%. However, our method of using 
smaller, source-free regions near the CDG should not be 
as sensitive as the plane-fitting method to residual large- 
scale gradients that may remain due to the dark current 
problems discussed in Section |3.1.1| and that may bias 
the derived sky level in the region of the CDG. 

We also use the difference between sky modes in the 
multiple regions to estimate the systematic e r ror in our 
sky estimate. As in McNamara & O'Connell (19921, we 
use the maximum difference between adopted sky modes 
(AS'inax = S'max " 5*111111) to define the systematic error 
in sky levels as AS". 



sys 



^ See http://cow.physics.wisc.edu/~craigm/idl/fitting/ 



AS'inax/2. When computing 
total errors in colors, systematic errors were added in 
quadrature with the statistical errors. 
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3.1.3. Image stacking 

After sky subtraction and before stacking, the pixel 
values were changed to units of electrons per second by 
multiplying by the gain and dividing by the exposure 
time. Next, the individual frames in each band were 
scaled to remove differences in intensity due to changes in 
airmass between exposures. The IRAF task mscimatch 
was used to perform this scaling. Mscimatch uses the 
measured intensities of unsaturated stars to determine 
the scalings; typically, 50-200 stars per frame were used. 
Since the frames for a given object were always taken 
within a short time of one another, scalings were typi- 
CE^Uy ^ 1%- We note that mscimatch was not used to 
determine zero-point offsets, as this step was done during 
sky subtraction. 

In order to stack the frames into a final master im- 
age, the frames must be tangent-plane projected onto 
the same pixel grid. To minimize artifacts due to projec- 
tion, cosmic rays were identified and removed using the 
craverage and fixpix tasks in IRAF. Mscimage was used 
to perform the projection using a sine interpolant. Next, 
in the rare cases where there were large point spread 
function (PSF) variations between frames, PSF match- 
ing was done to remove these variations. Finally, the 
frames were stacked using a median filter (or an average 
when only 2 frames were used) and a a-clipping routine 
with the mscstack task in IRAF. 

3.1.4. Calibration 

For the purposes of this paper, in which we com- 
pare colors measured between different radii of the same 
source, photometric calibration is unnecessary and was 
not done. While it is possible that we may miss signif- 
icant active star formation if it is distributed smoothly 
across the galaxy (resulting in possibly negative color 
gradients) , calibrated colors for similar systems (e.g., Mc- 
Namara fc 0'Connell|1992| are consistent with there be- 
ing little or no active star formation. Therefore, we con- 
sider it unlikely that significant star formation is present 
in any of the systems in our sample with negative color 
gradients. 

Usually, because our sources are scattered across the 
sky, corrections must be made to remove extinction due 
to dust in our own galaxy, which varies as a function of 
position on the sky. However, we are measuring radial 
gradients and hence comparing colors between two radii 
in the same source. Since the Galactic extinction does 
not vary significantly over the typical angular scale of 
our galaxies (r < 2 arcmin), no correction for Galactic 
extinction is necessary. 

3.2. X-ray Data 

AH systems were observed with the Chandra ACIS de- 
tector in imaging mode and the data were obtained from 
the Chandra Data Archive. Details of the observations 
are given in Table |2] 

The Chandra data were reprocessed with CIAO 3.3 
using CALDB 3.2.0 and were corrected for known time- 
dependent gain and charge transfer inefficiency problems. 
Blank-sky background files, normalized to the count rate 
of the source image in the 10 — 12 keV band, were used 
for background subtraction.* 



4. DATA ANALYSIS AND RESULTS 
4.1. Derivation of ICM Properties 

X-ray spectra were extracted in elliptical annuli with 
~ 3000 counts and were centered on the centroid of the 
cluster emission with eccentricity and position angle set 
to the average values of the cluster isophotes. In systems 
with clear cooling cusps, the X-ray centroid was fixed to 
the position of the core, determined by fitting a gaussian 
to the surface brightness distribution in that region. In 
systems without cooling cores, the X-ray centroid was 
determined by fitting an ellipse to the image at an in- 
termediate radius, usually at radii beyond any substruc- 
ture. However, in some of these systems (e.g., A520), 
due to the diffuse nature of the cluster emission or to 
substructure, it was very difficult to identify the cluster 
center. Therefore, the centroid positions for these sys- 
tems (noted in Table [l]) likely have large uncertainties. 
Spectra were extractea using dmextract in CIAO, and 
weighted response files were made using the CIAO tools 
mkwarf and mkacisrmf or mkrmf (mkacisrmf was used 
for all observations taken at the — 120 C focal plane tem- 
perature; mkrmf was used for all other observations). 

Gas temperatures and densities were found by depro- 
jecting the spectra with a single-temperature plasma 
model (MEKAL) with a foreground absorption model 
(WABS) using the PROJCT mixing model in XSPEC 
11.3.2, between the energies of 0.5 keV and 7.0 keV. 
The redshift was fixed to the value given in Table [T] and 
the foreground h ydrogen column density w as fixed to the 
Galactic value of Dickey & Lockman ( 1990 1, except in the 
cases of 2A 0335+096 and A478, when a significantly dif- 
ferent value was preferred by the fit. In these two cases, 
the column density in each annulus was allowed to vary. 
The density was then calculated from the normalization 
of the MEKAL component, assuming Tie = 1.2nH (for a 
fully ionized gas with hydrogen and helium mass frac- 
tions of X = 0.7 and Y = 0.28), as: 



'1.2 X 10i4(47rL>2) 



(l + z)2F 



(1) 



where rif, has units of cm^'', the luminosity distance (-D^) 
has units of cm, and the volume of the shell {V) has units 
of cm"^ . 

We derived the cooling times using the deprojected 
densities and tempe ratures found above and the cooling 
curves of Smith et al.^ (2001) , calculated using APEC^ 
(the resu lts do not ciiange significan tly if we use the 
curves of Bohringer & Hensler 19891. The pressure in 
each annulus was calculated as p = nkT, where we have 
assumed an ideal gas and n » 2ne- Lastly, the entropy 



is defined as in |Lloyd-Davies et al. (2000j) 



S = kTn-^/^, 



(2) 



and has units of keV cm^ . 

A primary objective of this study is to determine how 
the ICM properties relate to the optical colors. In par- 
ticular, we are interested in properties at the core of the 
cluster, where the cooling times are the shortest. How- 
ever, the systems in our sample span a wide range of 
redshift and core surface brightness, resulting in a wide 



See http: //asc .harvard. edu/contrib/maxim/acisbg/ 



® See http://cxc.harvard.edu/atomdb/ 
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range in the physical size of the innermost region (which 
we required to contain ^ 3000 counts). 

The situation can be improved somewhat using the 
fact that the emissivity and hence surface brightness for 
a fully ionized plasma is a strong function of the density 
(e ~ rig) to extrapolate the density closer to the core than 
is generally possible using spectral deprojection. We used 
this method for 15 of the 46 objects with optical data; 
for the remaining objects, this method did not result 
in a significant decrease in the central region over that 
used in spectral deprojection (due either to reaching the 
resolution limit of r « 1 arcsec or to substructure in the 
cluster emission that resulted in unphysical emissivities 
when the central region became too small). 
We used the onio n-peel deprojection method (e.g. 



Ameglio et al. 20071 to transform from X-ray surface 
brightness to density. We extracted the surface bright- 
ness profile between 0.5-7 keV in circular annuli with 
« 200-500 counts each, after background subtraction. 
When deprojecting the surface brightness profile, we ne- 
glected any emission from the cluster beyond the outer- 
most radius of the surface brightness profile, which can 
result in errors in the outer parts of the density profile if 
there is still significant cluster emission beyond this ra- 
dius. However, as we are interested only in the innermost 
portions of the profiles, this effect should not affect the 
central density that we derive. 

We then normalized the resulting density profile to 
the coarser density profile derived using the PROJCT 
model, and we used this normalized density profile to es- 
timate the central density. Errors on the derived central 
densities were estimated using a Monte Carlo technique, 
wherein the deprojection was repeated 100 times, with 
the surface brightness in each annulus drawn randomly 
from a poisson distribution given by the total counts in 
each annulus (before background subtraction). We also 
included errors in the central density due to uncertainties 
in the normalization by assuming, conservatively, that 
this error is equal to the error of the inner density derived 
from spectral deprojection. Finally, when calculating the 
central cooling time and entropy using this density, we 
assumed that the temperature and abundance of the gas 
are constant in the inner region used in spectral deprojec- 
tion and equal to the derived emission- weighted values. 
Table [3] lists the resulting central properties (and l-tr er- 
rors), derived as close as possible to the core (after the 
exclusion of any non-thermal point sources) . 

Lastly, to make comparisons of the ICM properties be- 
tween objects at a single physical radius, we used linear 
interpolation of the profiles in log — log space, using the 
mean radius of each annulus. We chose a radius of 12 kpc, 
as this radius was the smallest physical radius that could 
be achieved across most of our sample (excluding only 
A520; the mean radius of A1835 Hes at 12.7 kpc, but the 
properties at this radius should differ only slightly from 
those at 12 kpc). We also interpolated the cooling-time 
profiles in the same way to find the cooling time at the 
radius corresponding to the distance between the X-ray 
core and the CDG's core (denoted Ar and listed in Table 
[1]). Table [4] lists the interpolated ICM properties. 

4.2. Analysis of Optical Data 

We identified the CDG as the most extended galaxy 
within the field, which generally lies near the X-ray core. 



However, when the cluster contains more than one very 
large galaxy and it is unclear which is the dominant 
galaxy, we chose as the CDG the galaxy with the largest 
integrated if -band magnitude from the 2MASS catalog. 
Below we note any unusual characteristics of the CDG 
or image. 



A85: The CDG has extended optic al line emission ( Hu 
eFaTlfiJSSl [Fisher et al. II 1995[) 



3C 28: An excess of U V tiux has been detected in the 



CDG by I Wills et al | ( |2002| ) and was ascribed by 
them to active star formation. 

A262: The CDG has a large, central dust lane and a 
bright star nearby which were masked. 

Perseus: The Perseus CDG (NGC 1275) is well known 
for its blue, emission-line filaments, dust lanes, and 
the foreground high-velocity system, all of which 
were masked before anal ysis. Significant star for- 
mation was detected by McNamara fc 0'Connell| 
( |1989| ); |Romanishin| ( [l987P ; |Smith et al.| ( |1992p . 

2A 0335+096: There is a nearby bright star in the im- 
age that results in large areas of the galaxy being 
masked. Significant star forma tion was detected by 



Romanishin fc Hintzen (19881. 
A478: Cardiel et al. (1998) found spectral evidence for 



significant star formation in the CDG. 
A520: The cl uster appears to h ave no well-defined opti- 
cal core (Dahle et al. 12002), and instead has several 



concentrations of gal axies, none of wh i ch co incide 
with the X-ray core. Crawford et al. (1999) note 



that there are three dominant galaxies and chose 
the SW one as the CDG, as did we. 
M87: The optical jet and central region, associated with 

the AGN, were masked before analysis. 
HCG 62: The dominant galaxy of this compact group 
has a large, nearby companion galaxy which was 
masked before analysis. 
A1795: The CDG is well known to harbor a tail of blue 
emission that is thought to be due to star forma- 
tion triggered by compression of the gas by the 



radio source dMcNamara et al.| 1996a|b| |Pink ney] 

et al.|1996||Q'Dea~et al.|2004| . The R- and 7-b;md:^ 
images were att'ected by scattered light, which was 
masked as thoroughly as possible. 
A1991: As with A1795, the long-wavelength images 
were affected by scattered light, requiring masking, 
and th e seeing was poor. (McNamara fc O'Connell] 
(1989 ) report spectral evidence for a small SFR, 
and McNamara & O'Connell (19921 report a posi- 
tive U — 1 color gradient in the core. 
RXC J1504.1-0248: The [/-band image shows an elon- 
gated blue region in the core, extending 6 arcsec in 
length. 

A2065: This cluster has two dominant galaxies. We 
have analyzed the southern galaxy, whi ch appears 



to be asso ciated with the X-ray core (Chatzikos 
et al.''2006). 



RX J 1532.5+3021: Th e CDG is known to be very blue 
(Dahle et al. 2002). The U — I colormap shows a 
very blue region offset slightly to the south of the 
core and extending 4-8 arcsec in radius. 

Hercules A: The CDG has extended optical line emis- 
sion ( )Tadhunter et al.„1993j ). The CDG core is the 
southeastern of the two surface-brightness peaks. 
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TABLE 3 
Central ICM Properties. 







lie 


kT 


Entropy 


toool 


System 


(arcscc) (kpc) 


(cm-^) 


(kcV) 


(keV cm^) 


(10« yr) 



A85 


2.4 


2.6 


loT+o oof 


2.06l°o:f7 


9 1+0** 

■ — . 9 


1 7+0.5 
1- -0.5 


3C 28 


2.2 


7.2 


053+""^'^ 

U.UOO_Q QQQ 


1.29+°o;^ 


9.1+?? 


2.8i;-° 


A133 


3.0 


3.3 


O.U48_o 


1 77+0.09 
' ' -0.06 


13.5+;i 


9 9+0.9 

^■■^-0.6 


A223 


2.0 


6.7 


U.UJ-y_Q 003 


5.50«:«° 




24.0l^5;"2 


A262 


2.0 


0.7 


067+° "^'^ 

u.uu I _o.019 


0.86l°;°o; 


5.2tl° 


n s+O-2 
0.8_o.2 


A383 


1.5 


4.6 


116+0010 

U.J-J-U_o 012 


Z.UO_Q 26 




1 s+o-'' 

l-°-0 6 


AWM 7 


2.2 


0.8 


140+0 015 

U.l'iU_o 013 


. ..0+0. 07 
-'-■-'-^-0.07 






Perseus 


11.8 


4.2 


1 i^n+ooos 

U.10U_o 005 




15.5t;- 


9 c: + 0.4 


2A 0335+096 


3.4 


2.4 


100+°"^^ 
u.iuu_o oil 


. 09+0.07 
-'-■OZ_Q Qg 


6 1+°° 


1 4 + 0-6 


A478 


1.5 


2.6 


197+0 014 
u.iy/ -0.015 


2-68«;28 


7 g + 0-9 
' ■ y — . 9 


-1 + 0.4 
-'--'-'-0.3 


A496 


1.0 


0.6 


198+°"^^ 
u. j-yo_Q 025 


1 20+° 1^ 
-l--^U_o 15 


r; + 0.5 
''■°-0.5 


0.7l°o;^3 


A520 


7.7 


25.2 


005+0 001 
U.UUO_o 001 


7 07+^'2'' 

' ■'J' -0.90 


253.8if3°2* 


118-0;°8.o 


MS 0735.6+7421 


2.9 


10.1 


067+°"''^ 

u.uu / _o.002 


10+0.22 
■'■l^-O 24 




4 7+0-7 
^- ' -0.6 


PKS 0745-191 


2.5 


4.8 


143+0.010 

^■-"-^■^-O-OIO 


2 62+°-^^ 

'^■'^'^-0.36 




9 1+0.6 
^--1-0.5 


Hydra A 


2.2 


2.4 


149+0 01-1 

^■-^^^ — 0.019 


2.59t°o:S 






Zw 3146 


1.6 


6.8 


n 1 77+0.007 
^■-L ' ' -0.007 


3.09t°- 


9 8+°° 


1 9 + 0-3 
-^-0.3 


A1068 


2.0 


4.9 


149+0 006 
u.r'iy-o.oii 


2 16+°''2 

'^■10-0.09 


7 7+1.5 
' ■ ' -0.5 


i-4lS:=3 


A1361 


2.0 


4.2 


044+°"^° 

U.U'±'±_o 010 


2 5g+° io 

^.oy_0 22 


2o.7j:;^;l 


5.7i;-; 


A1413 


3.8 


9.5 


039+° ™^ 
u.uoy_o 001 


6 14+135 
"■-L^-1.00 


53.4^11/ 


,0 1+2.9 
J^J-J^-2.6 


M87 


5.9 


0.5 


191 + 0.009 
^■-^^-^ — 0.009 


94+0 02 
^■y^-0.02 




4+°-^ 


HCG 62 


2.0 


0.6 


072+o oi'i 


n C7+0.01 

0-671o.oi 


3.9t°o:°6 


6+°-^ 


A1650 


3.0 


4.7 


043+0 002 

U.U^O„o.o02 


^■^O_o 90 


34.41 1^ 


-1+2.8 
8.1-1 8 


Coma 


9.8 


4.6 


n nnH+0 003 
0.008_o.oo3 


_ +3.45 
^- '^-1.61 


66.9lf2"6 


09 0+^3.8 
OZ.U_3i 3 


A1795 


2.0 


2.4 


069+° °^° 

u.uu»_0 Qig 


2 74+0.57 

^- '^-0 37 






A1835 


3.2 


12.7 


110+0 003 
U.J-J-U_o Q03 


4 03+°-" 

'±.uo_o 27 




3-Ot°o1 


A1991 


2.2 


2.5 


077+0 011 

U.Ui I -0.011 


Q 77+0.02 

^- ' ' -0.02 


4.2l°o:l 


0-6to.l 


MS 1455.0+2232 


1.3 


5.3 


095+° °^'^ 
u.uyo_o 020 


4 40+^ ='° 

^■^U-0.90 


21.1^5° 




RXC ,11504.1-0248 


1.1 


3.8 


n isn+o oii 

0.i8U_o 012 


6.87l^;«^ 


91 r;+12.0 




A2029 


0.6 


0.9 


n 070+0. 036 
U.OiO_o 031 


o + - 33 
^■88-0 22 


5.6l°o:I 




A2052 


3.9 


2.7 


0-017t°o:™02 


7i+°-04 

U. ' J-_o 08 


10.511° 


r + 7.a 
•^-^-1.8 


A2065 


2.0 


2.7 


o.o37«:S; 


1 97+0.35 
— 0.21 


17.S+_ti 




RX J1532. 8+3021 


2.0 


9.6 


107+0 009 

u. J-U I -0.009 


09+0.29 
''■'^'^-0.28 


14.8;;° 


+ 0.6 
•^-"-0.5 


A2218 


3.9 


11.7 


n nnR+0 003 
0.008_o.oo3 


3 44+°-00 

•^■^^-0.56 


00 0+32.1 

00.0_25 9 


■^4 r + 22.2 
o'i-0_18.1 


Hercules A 


2.0 


5.3 


n n7«+o oio 
0.078_o.oio 


2 04+°-i^ 

^•"^-0.21 


11 2+1'^ 




A2244 


2.0 


3.5 


046+° ™*' 

U.U^D_0.008 


4 58+°-°2 

1-Oo_o 00 


35. 6^^ J 


a K + 3-1 
6-6-2.1 


NGC 6338 


1.5 


0.8 


236+0 022 
U.^JO-0.018 


99+° o^ 

•^■^^-0.06 


2.6;°„:i 


5+01 


RX J1720. 2+2637 


1.7 


4.8 


099+° °^^ 

U.Uy»_0 Oig 


9 n7+0-83 
•^■^ ' -0.61 


13.81*3- 




MACS J1720. 2+3536 


2.2 


11.7 


076+° ™^ 

U.U(D_o 006 


cc + 0-34 
3-66-0.25 


20.4+ig 


r + 1-0 
•3-0-0.8 


A2261 


2.0 


7.1 


043+° °^° 
u.uio_o 010 


5 7fi+l-35 
^- '"-0.92 


46 6+^^ * 


7 + 4.3 
'-■3-3.0 


A2319 


4.4 


4.8 


n n':>»+o oi4 
0.02S_o 014 


11 32+''-*'' 
li-J^-3.41 


,99 9 + 90.4 
±z.i.z_55 Q 


24 4+^''-° 


A2390 


0.6 


2.3 


199+0 010 
u. j-yy_o 015 


3 'iq+1-35 
J-'3»-0.86 


io.o;*2.°6 


2 0+°° 


A2409 


3.0 


7.6 


013+0 008 
U.U10_o 008 


4 82+°-°** 

*-''^-0.68 


88.5tfo1 


21 4+1**-° 

^l-*-16.0 


A2597 


3.4 


5.5 


075+0 005 
U.U(0_o 005 


1 c^ + 0-12 
l-66_o.l4 


Q 0+0.8 

a-J-o.o 


2 6+°° 


A2626 


2.0 


2.1 


063+0 01'l 
U.UD.3_0.014 


2 54+°-'='^ 

^•'^^-0.64 


16 0+^ '' 


4 1+2-1 

— 2 2 


A2657 


3.0 


2.3 


017+0 000 

U.Ul 1-0.009 


04+1.39 
3.84_o.87 


58 2+20-'' 
00-^-24.6 


i8.4i;t? 


A2670 


3.0 


4.3 


026+° °^^ 

U.UZD_o.ol2 


00 + 0.90 
J-'JO_o.65 


00 + 17.0 
38.d_i5 9 


6-6i^:'o 



^ Mean radius of the inner-most region; when elliptical regions -were used, the equivalent redius 
(r — [afc]i/2^ where a and h arc the semi-axes) is given. 



KTIAA: A nearby bright star required large areas of 
masking. 

RX J1720.2+2637: The CDG has a blue central re- 
gion some 13 arcsec in radius. 

MACS J1720.2+3536: The U ~ 1 colormap shows 
blue emission extending to a radius of 2.5 arcsec. 

A2390: The [/-band image shows a very blue, elongated 
region in the core, extending SE to NW some 4-5 
arcsec in length. 

A2 5 9 7: The CDG possesses knots of star fo rmation 
dKoekemoer et al.||2002| [O'Dea et aL]|2004 l. 



4.2.1. Surface Brightness Profiles 

The master frames, created by stacking the individual 
frames for each object (see Section 3.1.31, were first reg- 
istered to a common coordinate grid using wregister in 
IRAF. Next, the master frames were PSF-matched using 
the psfmatch task to the lowest-resolution frame of the 
set (typically the {/-band frame). This step is necessary 
before spatial comparisons between frames taken in dif- 
ferent bands can be made. However, since the bands are 
fairly broad {FWHM - 1500 A) and the stars used for 
the matching may have different colors than the CDG, 
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TABLE 4 
Interpolated ICM Properties. 
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The Ar radius is the projected distance between the X-ray core and the CDG core (see 
Table [TJ. 

^ The Coma data lack sufficient counts to derive a cooling time at the CDG's radius. 



some effect due to PSF mismatches -will remain in the 
data. 

Surface brightness profiles -were constructed from the 
PSF-matched, registered master frames using the el- 
lipse task in IRAF. This task -works by fitting ellipses 
to isophotes of the surface brightness distribution of a 
galaxy and calct ilating the mean sur face brightness along 
the ellipse (see Jedrzejewski 1987). Clipping routines 
are used to reject stars, superimposed galaxies, or other 
source of contamination. Additionally, such contaminat- 



ing objects as -were visible -were masked out by eye before 
fitting. The ellipticities and position angles -were allo-wed 
to vary from ellipse to ellipse to reflect changes in the 
underlying galaxy. Typically, the surface brightness dis- 
tributions of CDGs bec ome more elliptical at larger radii 
(e.g., ,Patei et al. 20061. The eUipse centroids were not 
allowed to vary, to make comparisons between profiles 
in different bands possible. The centroids were deter- 
mined from the long- wavelength images {R or /), which 
are generally smooth and relaxed (reflecting the old stel- 
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Fig. 1. — U — R {left) and U — I {right) color versus the equivalent radius (r = [afe]^/^, where a and b are the semiaxes), ordered by 
decreasing gradient. The colors for each object have been shifted by arbitrary values. The error bars shown for each point are the statistical 
errors; the dashed lines show the total (statistical plus systematical) errors. The best-fit gradients (see Section [4.3[ l are overplotted as solid 
lines between the inner and outer radii used in the fit. 

lar populations of the galaxy) . Centroid positions were 
calculated using imexamine in IRAF, which finds the 
centroid using Gaussian fits to the radial surface bright- 
ness profile. 

The ellipse centers, ellipticities, and position angles 
were fixed to those resulting from a fit to the R- or I- 
band image, whichever was of higher signal-to-noise, to 
make it possible to compare profiles in different bands. 
We note that fixing the [/-band ellipse properties to those 
of a longer-wavelength band will tend to reduce the U- 
band surface brightness along the ellipse if the [/-band 
emission has significantly different preferred ellipticities 
or position angles. An example of such a case is A2390, 
which has a bar of very blue emission near its center that 
is clearly different in shape from the J-band emission in 
the same region. The net effect of our procedure is to 
dilute somewhat the signatures of star formation in our 
profiles; however, we estimate that the effect is generally 
small. The resulting color profiles for the optical sample 
are shown in Figure [l] 

4.2.2. Comparison with Other Studies and between 
Observing Runs 

As a check for systematic errors in our reduction and 
analysis, we can compare our color profiles both to pro- 
files from the literature and between different observing 
runs. In Figure [2j we compare our U ~ I profiles with 




Fig. 2 . — Comparison of our U — I co lor profiles {triangles) to 
those of |McNamara &: 0'Connell| | [T992[ | {circles) for the objects 
common to both samples. 



those from Figure 5 of McNamara & O'Connell (19921. 
In general, the profiles agree well within the errors over 
the inner 10-20 arcsec. At the extreme inner radii there 
is often some discrepancy, probably due to seeing effects 
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Fig. 3. — Comparison of U — R color profiles {left) and U — I 
color profiles (right) from different observing runs: A2390 and RX 
J1720.2+2637-MDM (triangles) and KPNO (circles), A2261 and 
A2244 - MDM September 2005 (triangles) and MDM May 2006 
(circles). 



at a 5, 2 arcsec, as the images of McNamara & O'Connell 
( 1992 1 were not matched between bands for different see- 
ing. At the extreme outer radii (a > 10 — 20 arcsec) 
our profiles do not match well, as those of Mc Namara] 
& O'Connell ( |1992 1 often show a steep positive gradient 
that ours do not (particularly apparent in the profiles of 
A1991 and A262). This effect may be due to mise sti- 
mation of the sky by McNamara & O'Connell (19921, as 
they used a much smaller CCU (I'OV^ 2.5 x 2.b arcmin 
compared to our smallest FOV of 9.4 x 9.4 arcmin) which 
may have resulted in sky boxes with some contamination 
from the central galaxy's halo. Otherwise, the only sig- 
nificant discrepancy is in the 3 < a < 6 arcsec region of 
the profile of A1795. We believe the bluer colors that 
McNamara fc O'ConneU] ([1992) find in this region are 
due to ditterences in the aiinuii used to measure the sur- 
face brightness, as the J7-band image shows a linear blue 
feature in this region, the color of which would be diluted 
if fairly circular annuli are used (as we have done). 

As a further check of consistency, we also compare pro- 
files taken from observations made during different runs 
and at different telescopes. In Figure |3] we show profiles 
for three objects that we observed twice. A2390 was ob- 
served with the MDM 2.4-m telescope and the KPNO 4- 
m telescope; A2261 and A2244 were observed twice with 
MDM. The profiles agree well, with the exception of the 
inner 2 arcsec of the A2390 profiles, which is probably 
due to seeing differences between the two observations. 

4.3. Color Gradients and A Colors 

Color gradients were derived for the U — I and U — R 
color profiles using a least-squares fitting routine in the 
MPFIT package. The color gradient is defined as the 
change in color (in magnitudes) over the corresponding 




-0.5 0.0 0.5 1.0 1.5 
G(U-R) (mag/dex) 

Fig. 4. — The U — I color gradient versus the U—R color gradient. 
The dotted line shows equality between the two gradients, while 
the solid line shows the best-fit straight line. 

change in log(r) (in dex); e.g., for U — I: 



G{U - I) 



d{U - I) 
rflog(r) ' 



(3) 



The following function was fit to the data in color-log(r) 
space: 

U-I = G{U - I) log(r) + 6, (4) 

where G{U — I) is the color gradient and h is the inter- 
cept. The data were weighted by their total (systematic 
plus statistical) error, and errors on the gradient were 
returned by the covariance matrix of the best fit. To 
eliminate possible effects from the PSF on the fitted gra- 
dient, the fits were restricted to radii greater than twice 
the radius of the FWHM of the PSF. 

The outer radius used in the fit was set differently 
for objects with blue cores and for those without. For 
blue objects, the outer radius was set to the approxi- 
mate point at which the excess blue emission ends, so 
that the resulting gradient would trace the star-forming 
region only. The edge of the blue emission was deter- 
mined by examining the U — I or U — R colormaps of 
the galaxy, made by dividing the [/-band image by either 
the R- or /-band image. For red objects, the outer ra- 
dius was set to the radius at which the total errors reach 
0.5 mag to avoid regions where sky-subtraction errors 
significantly affect the profile (however, since the data 
were weighted by their errors, this choice has little effect 
on the measured gradients). See Figure [T|for plots of the 
color profiles with the best-fit gradients overlaid. Table|5] 
lists the color gradients and the radial range over which 
they were derived. 

The U — R and U — I color gradients are very similar, 
with the U — I color gradients being on average slightly 
more positive (due to a lower contribution from the star 
formation to the flux in the / band). To illustrate this, 
we plot in Figure [4] the two gradients against one an- 
other. Since the errors in the two gradients are correlated 
through the errors in the U -band surface-brightness pro- 
files, we perform a li near regression using t h e BCE S least- 
squares method of Akritas & Bershady (19961, which 
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TABLE 5 

CDG Color Gradients and A Colors. 





Color Gradients 






A Colors 






GiU - R) 


GiU - 1) 


Range 


A{U - R) 


A(U - I) 


Range 


s 

ye, cm 


(mag dcx~ "'" ) 


(mag dex "'" ) 




(mag) 


(mag) 


(kpe) 


AS 5 


0.02 ± 0.04 


0.06 ± 0.05 


1.6-9.5 


-0.03 ± 0-05 


0.00 ± 0.05 


5-10 










-0-10 ± 0-10 


-0.08 ± 0.10 


5-20 


3C 28 


0.72 ± 0.26 


0.70 ± 0.30 


1.3-5.1 


0.26 ± 0.14 


0.24 ± 0.15 


5-10 


A133 




-0.21 ± 0.04 


1.8-39.6 




-0.09 ± 0.05 


5-10 












-0.19 ± 0.10 


5-20 


A223 


-0.34 ± 0.31 




1.5-8.9 


-0-14 ± 0-18 




5-10 










-0-06 ± 0-39 




5-20 


A262 




-0.32 ± 0.03 


2.0-176.4 




-0.06 ± 0.10 


5-10 












-0.19 ± 0.25 


5-20 


A383 


0. 15 ± 0.28 


0.10 ± 0.25 


1.8-5.3 








AWM 7 




-0.18 ± 0.03 


1.1-84.8 




-0.09 ± 0.18 


5-10 












-0.17 ± 0.45 


5-20 


Perseus 




1. 19 ± 0.07 


6.5-47.0 




0.38 ± 0.05 


5-10 












0.59 ± 0.17 


5-20 


2A 0335 + 096 


0.14 ± 0.08 


0.25 ± 0.09 


1.5-4.2 


-0.12 ± 0.08 


-0.22 ± 0.08 


5-10 


A478 


0.84 ± 0.25 




1.0-5.1 


0.16 ± 0.53 




5-10 


A496 


-0.15 ± 0.05 


-0.08 ± 0.06 


1.7-18.7 


-0.05 ± 0.08 


-0.04 ± 0.09 


5-10 










-0.08 ± 0.22 


-0.11 ± 0.24 


5-20 


A520 




— 0.20 ± 0.28 


1.2-9.0 




-0.13 ± 0.16 


5-10 












-0.13 ± 0.48 


5-20 


MS 0735.6+7421 


— 0. 17 ± 0.24 


-0.08 ± 0.28 


1.0-4.9 


-0.03 ± 0.24 


0.00 ± 0.24 


5-10 


PKS 0745-191^ 




1.496 ± 0.13 


2.0-6.0 








Hydra A^ 




1.0 ± 0.16 


1.0-5.0 








Zw 3146*^ 


3.4 ± 0.18 




2.1-8.0 








Aloes'^ 


1.7 ± 0.18 




1.5-4.0 








A1361 


-0.15 ± 0.06 


-0.07 ± 0.06 


2.0-13.1 


-0.04 ± 0.03 


-0.03 ± 0.03 


5-10 










-0.06 ± 0.06 


-0.00 ± 0.07 


5-20 


A1413 


-0.29 ± 0.15 


-0.22 ± 0.21 


1.7-9.6 


-0.05 ± 0.09 


-0.03 ± 0.10 


5-10 










-0.23 ± 0.18 


-0.19 ± 0.21 


5-20 


M87 


-0.27 ± 0.08 


-0.20 ± 0.07 


11.5-96.6 


-0.15 ± 0.43 


-0.16 ± 0.42 


5-10 


HCG 62 


-0.33 ± 0.04 


-0.24 ± 0.05 


2.3-42.1 


-0.07 ± 0.11 


-0.14 ± 0.13 


5-10 


A1650 


-0.14 ± 0.11 


-0.07 ± 0.14 


2.1-19.7 


-0.05 ± 0.11 


-0.02 ± 0.11 


5-10 










-0.20 ± 0.17 


-0.14 ± 0.20 


5-20 


Coma 


-0.14 ± 0.02 




1.7-94.4 


-0.04 ± 0.08 




5-10 










-0.06 ± 0.18 




5-20 


A1795 


0.50 ± 0.10 


0.68 ± 0.11 


2.1-9.1 


0.05 ± 0.08 


0.14 ± 0.08 


5-10 










0.07 ± 0.18 


0.14 ± 0.20 


5-20 


A1835^ 


4.5 ± 0.18 




2.0-7.0 








A1991 


-0.35 ± 0.06 


-0.23 ± 0.08 


2.6-16.9 


-0.11 ± 0.05 


-0.08 ± 0.06 


5-10 










-0.13 ± 0.10 


-0.02 ± 0.12 


5-20 


MS 1455.0+2232 


-0.06 ± 0.23 


0.02 ± 0.25 


1.5-5.8 








RXC J1504. 1-0248 


0.94 ± 0.17 


1.03 ± 0.16 


1.7-5.8 






5-10 


A2029^ 




-0.269 ± 0.014 


3.0-10.0 








A2052^ 




-0.069 ± 0.021 


5.0-35.0 








A2065 


-0.38 ± 0.08 


-0.31 ± 0.08 


1.4-29.2 


0.01 ± 0.10 


0.04 ± 0.10 


5-10 










-0.20 ± 0.18 


-0.11 ± 0.18 


5-20 


RX J1532. 8 + 3021 


1.04 ± 0.40 


1.56 ± 0.53 


2.0-5.4 








A2218 


-0.41 ± 0.11 


-0.40 ± 0.12 


1.4-8.7 


0.02 ± 0.10 


0.02 ± 0.10 


5-10 










-0.24 ± 0.11 


-0.25 ± 0.11 


5-20 


Hercules A 


0.34 ± 0.25 




1.6-8.8 


0.13 ± 0.11 




5-10 










0.20 ± 0.32 




5-20 


A2244 


-0.20 ± 0.12 


-0.02 ± 0.13 


1.9-18.1 


-0.12 ± 0.09 


-0.05 ± 0.10 


5-10 










-0.28 ± 0.17 


-0.17 ± 0.17 


5-20 


NGC 6338 




-0.12 ± 0.01 


1.6-102.2 




-0.04 ± 0.03 


5-10 












-0.14 ± 0.09 


5-20 


RX .11720.2 + 2637 


0.28 ± 0.14 


0.46 ± 0.17 


1.5-11.7 


0.14 ± 0.10 


0.16 ± 0.11 


5-10 










0.27 ± 0.18 


0.41 ± 0.19 


5-20 


MACS ,11720.2 + 3536 


.84 ± .85 


1 .06 ± .88 


1.6-2.2 








A2261 


-0.42 ± 0.19 


-0.35 ± 0.19 


1.4-7.4 


-0.10 ± 0.10 


-0.10 ± 0.10 


5-10 










— 0.31 i 0.17 


— 0.26 i 0.17 




A2319 




-0.13 ± 0.05 


1.3-18.0 




— 0.06 ± 0.07 


5-10 












-0.14 ± 0.13 


5-20 


A2390 


1.35 ± 0.67 


1.00 ± 0.76 


1.9-4.0 


0.77 ± 0.17 


0.67 ± 0.17 


5-10 










1.41 ± 0.54 


1.39 ± 0.56 


5-20 


A2409 


— 0.34 ± 0.26 


-0.16 ± 0.28 


1.5-6.3 


-0.12 ± 0.17 


-0.08 ± 0.17 


5-10 


A2597 


0.55 ± 0.14 


0.72 ± 0.14 


2.1-8.6 


0.23 ± 0.10 


0.31 ± 0.10 


5-10 










0.23 ± 0.32 


0.29 ± 0.32 


5-20 


A2626 


-0.09 ± 0.03 


-0.02 ± 0.04 


1.1-18.8 


-0.01 ± 0.07 


-0.01 ± 0.07 


5-10 










-0.03 ± 0.13 


0.03 ± 0.13 


5-20 


A2657 


-0.14 ± 0.06 


-0.08 ± 0.07 


1.4-9.5 


-0.02 ± 0.15 


0.01 ± 0.15 


5-10 










-0.01 ± 0.39 


0.10 ± 0.39 


5-20 


A2670 




-0.26 ± 0.03 


1.4-32.2 




-0.08 ± 0.03 


5-10 












-0.18 ± 0.06 


5-20 


^ Gradient taken from 


McNamara & O'C 


onnellNl992| 










^ Gradient taken from 


McNamara| jl995| 


W^e have adopted the average 


G(U - I) error fo 


r our sample. 




Gradient taken from 


unpublislied data 


(B . R. McNamara, private communication) . We 


lavc adopted the 


average 


G(U - R) error for ou 














^ Gradient taken from 


McNamara &z O'C 


onnellNl992l. We have adopted the avcraee G(U 


- R) error for our 


sample. 


Gradient taken from 


McNamara et al.| 


2006(. We have adopted the av 


erage G (U - fi) e 


ror for our sampl 





properly accounts for correlated errors. The resulting fit 
(with l-cr errors) is: 

G{U-I) = (0.090±0.015) + (1.01±0.14)G'(C/-i?). (5) 

The U ~ I color gradients are therefore consistent with 
a simple offset from the U ~ R color gradients. We note 
that there are no objects for which the difference between 
the two gradients is large enough (> 2a) to indicate sig- 
nificant contamination from emission lines in the R or I 
images. 



The gradients discussed above do not by themselves 
give an indication of the radial extent of the star forma- 
tion. For two objects with the same gradient, one may 
have star formation over a much larger physical radius. 
Therefore, it is useful to derive the change in color be- 
tween two physical radii. Objects with star- forming re- 
gions that are physically smaller will tend to have smaller 
color changes than those objects with more extended 
star formation. To derive the color changes, denoted 
A([/ — i?) and A(J7 — /), we use simple linear interpo- 
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lation to estimate the colors at 5, 10, and 20 kpc. The 
A{U — R) color changes are then: 

A{U -R) = {U- RUter - (C/ - i?)i„„cr (6) 

Table |5] lists the A colors for the optical sample. 

The detectability of blue emission depends on sev- 
eral factors, such as the signal-to-noise ratio of the im- 
ages, the seeing or PSF size (as we exclude unresolved 
emission), and the redshift of the source. Since this 
PSF size limits our spatial resolution, an obvious bias 
is that we are sensitive to blue emission that extends 
over smaller physical radii in lower-redshift systems than 
in higher-redshift systems. This bias could be impor- 
tant if the red systems are at systematically higher red- 
shift than the blue systems, such that extended blue 
emission could be missed preferentially in the red sys- 
tems. However, the mean redshift of the blue systems 
(< z >— 0.159) is approximately twice that of the red 
systems (< z >— 0.073). The lower mean redshift of 
the red systems is likely due to X-ray sele ction effects, 
because, as we demonstrate in Section 5.1 the red sys- 
tems often lack short central cooling times and cool cores, 
making them more difficult to image in X-rays at a given 
redshift. 

5. DISCUSSION 

5.1. The Relation between Star Formation and Cooling 

Time 

The cooling time of the ICM is a critical measure of 
its thermal state. If the cooling time is sufficiently short 
(on the order of the age of the system or less), signif- 
icant condensation of gas should occur unless heating 
balances the cooling entirely. In this section, we investi- 
gate whether the cooling times we derived from Chandra 
data are related to the presence of excess blue emission 
in the CDG. 

To this end, in Figure |5] we plot the optical color gra- 
dients of the CDG against the central cooling times of 
the cluster's ICM. We plot the U — I color gradient when 
available; otherwise, we plot the U — R gradient, as they 
are almost equivalent (see Section 4.3). The left panel 
shows the emission-weighted central cooling time derived 
as close to the core as our data allow us to achieve (see 
Table l3l . Since the systems in our sample vary greatly 
in redsnift, we plot in the middle panel the cooling time 
at a single physical radius of 12 kpc (see Table |4]). At 12 
kpc, the very short cooling times (those below 2 x 10^ 
yr) measured at the center disappear, demonstrating that 
their short cooling times relative to the rest of the sample 
are due to resolution effects. 

For comparison, we also plot the A colors against the 
cooling times in Figure [6] In this figure, the overall re- 
lationship between cooling time and the presence of star 
formation is unchanged from that shown in Figure |5] If 
the star formation rate is governed by the local value of 
the cooling time at each point throughout the ICM, then 
we should expect the A colors to increase smoothly as 
the cooling time at 12 kpc decreases. The abrupt turn 
on of star formation at tcooi — 4-8 x 10* yr in Figure [b] 
does not support such behavior. 

It is apparent from Figure [5] that positive gradients, 
which are indicative of active star formation, occur only 
in objects with central (r < 12 kpc) cooling times be- 
low 7-8 X 10* yr. Such a direct relationship between 



star formation and the local state of the intracluster gas 
strongly implicates the gas as the cause of the star for- 
mation. The simplest possibility is that cooling gas fu- 
els the star formation directly. Alternatively, dense gas 
might have fueled an AGN outburst that triggered the 
star formation. However, it is also apparent that the 
presence of a short central cooling time does not guar- 
antee the presence of a positive color gradient, as several 
of the objects with cooling times below ~ 8 x 10* yr do 
not have blue cores (although the majority do). 

To emphasize that the systems with blue central col- 
ors are different on average from those with red colors, 
we show in Figure [7] histograms of the cooling times at 
12 kpc for the objects in our optical sample with cool- 
ing times below 4 x 10^ yr. We have divided the sample 
into those systems with positive color gradients (blue sys- 
tems) and those with negative gradients (red systems). 
We include MS 1455.0-^2232, which has a positive U-I 
gradient and a negative U — R gradient, in the blue 
sample. Compared with the red systems, the blue sys- 
tems clearly prefer shorter cooling times. A Kolmogorov- 
Smirnov (K-S) test gives a probability of just 1.3 x 10~^ 
that the red and blue samples are drawn randomly from 
the same parent distribution of cooling times. This re- 
sult is a strong indicator that the cooling flow and star 
formation in the CDG are connected. 

The large disparity between the mass of gas within the 
cooling regions of the cluster and the mass of cold gas 
and new stars, coupled with the sharply declining cooling 
times with decreasing radius, precludes the possibility 
that the central cooling time is significantly affected by 
the cold gas. Thus, the tight relationship between cooling 
times and color gradients make it implausible that the 
stars are formed from infalling cold gas. 

One possibility is that thermal instabilities in the hot, 
cooling gas are driving the star formation. Under the 
right conditions, the motion of a cooling blob of gas will 
become unstable, leading eventually to catastrophic cool- 
ing. These conditions are met roughly when the growth 
rate of the instabilities due to cooling exceeds the coun- 
teracting rate of any damping, such as that from viscous 
damping (which acts to slow the blob's radial oscilla- 
tions) or thermal conduction (which acts to reduce the 
blob's temperature contrast with its surroundings) . The 
stability condition in any region can be determined from 
the X-ray data, and a detailed treatment of this scenario 
will be presented in a forthcoming paper. 

5.2. Star Formation and Other ICM and Calaxy 
Properties 

We also investigated the dependence of the presence of 
excess central blue emission on the temperature, density, 
and entropy of the ICM. We plot the CDG color gradi- 
ents against each of these properties, derived at a single 
physical radius of 12 kpc, in Figure [8j The temperature 
shows no clear relation to the color gradient. However, it 
appears that high densities (> 0.03 cm^'^) are required 
for blue colors, in agreement with the cooling time find- 
ings discussed earlier (since tcooi "-cT^)- 

Additionally, we find that positive gradients occur only 
in objects in our sample whose entropies at 12 kpc are 
< 30 keV cm^ (see the right-hand panel of Figure 
The range in entropy levels seen in this plot for the blue 
objects is consistent with that found by [Donahue et al. 
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Fig. 5. — Color gradient versus the central cooling time {left panel), the cooling time at a radius of 12 kpc {middle panel), and the 
cooling time at the location of the CDG {right panel). U — I gradients {empty symbols) are used when available; U — R gradients {filled 
symbols) are used when U — I gradients are unavailable. The horizontal, dashed line marks the division between systems with positive 
color gradients (blue, star-forming systems) and those with negative color gradients (red systems). The vertical line in the middle panel 
marks tcool(12 kpc) = 8 X 10** yr. 
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Fig. 6.— A{U - R) and A{U - I) between 5-20 kpc versus the 
cooling time at 12 kpc. Symbols are the same as those in Figurejs] 



([2006 ) for a sample of 9 cooling- flow clusters (all o f which 



are shared with our sample) . iVoit & Donahue ( 2005 1 



show that such entropy levels of ~ 10-30 keV cm^ are 
consistent with those expected from the effects of AGN 
heating on a simple cooling entropy profile, implying an 
indirect connection between feedback and star formation. 

The great majority of CDGs host a radio source (?). 
Therefore, it is useful to examine whether the presence 
of a radio source (and hence an AGN) is related to the 
presence of star formation. To identify AGNs, we have 
searched available radio catalogs (e.g., the NVSS and 
FIRST catalogs) for radio sources associated with the 
CDGs in our sample. We find detections at 1400 MHz for 
42 of our 46 CDGs. The four sources that lack detections 
are A1413, A2218, AWM 7, and A2261, all of which have 
negative central color gradients (and consequently no ev- 
idence of recent star formation) and long central cooling 
times. However, the remaining 21 systems with negative 



Fig. 7. — Histogram of the cooling times at 12 kpc for objects 
with positive color gradients (blue systems, gray region) and neg- 
ative gradients (red systems, hatched region). Three red systems 
with cooling times greater than 4 X 10^ yr are not shown. 

gradients have radio sources and, therefore, the presence 
of an AGN does not appear to depend stongly on the 
presence of star formation in the CDG (and vice versa). 
A detailed study of the relation between the radio-source 
and star-formation properties is beyond the scope of this 
paper and will presented in a future paper. 

Lastly, we investigated whether the color gradient has 
some dependence on the total luminosity of the galaxy. 
We plot in Figure [9] the color gradient versus the total 
X-band luminosity of the CDG. Apparent, total if -band 
magnitudes were taken from the Two Micron All Sky 
Survey (2MASS) catalog.^ The apparent magnitudes 
were corrected for Galactic extinction with the values 



of Schlegel et al. (19981 and corrected for redshift {K- 



See |http: //www. Ipac . caltech.edu/2mass/ 
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Fig. 8. — Color gradient versus the temperature, density, and entropy at 12 kpc. Symbols and the dashed line are the same as those in 
Figure [5] 
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Fig. 9. — Color gradient versus the total if-band absolute mag- 
nitude of the galaxy. Symbols and the dashed line are the same as 
those in Figure [5] 



correc t ed) a nd evolution using the corrections of Pog- 
gianti (jl997|). Lastly, the magnitudes were converted to 
absolute magnitudes using our assumed cosmology and 
the redshifts listed in Table [l] (which also gives the re- 
sulting absolute magnitudes). It is clear from Figure |9] 
that the color gradient shows no clear relation to the to- 
tal luminosity of the galaxy; positive gradients are found 
across almost the full range of luminosities. Addition- 
ally, a clear dichotomy between quiescent red galaxies 
and blue star-forming galaxies with steep gradients is 
evident. 

5.3. Cooling, Star- Formation, and AGN-Heating 
Timescales 

In the simple AGN feedback model, the timescales over 
which cooling, star formation, and heating occur should 
be related. The presence of shocks, in pa rticular, sug- 
gest that AGN heating is intermittent (e.g., Nulsen et al 



2005a|b| |\yilson et al.|[2006l [Fabian et ai:][2006l [FOTSIaH] 
et al.|2007 1 . In such a model, the AGN must create cav- 
ities frequently enough to prevent large amounts of net 
cooling and to maintain the gas at the observed temper- 
atures. This condition can only be met when the average 
time between outbursts is less than the average central 
cooling time. Additionally, cooling must occur over a 
long enough period for significant cool gas to accumu- 
late and star formation to occur, after which the signa- 
tures of star formation may persist even though cooling 
has ceased. The probable interrelation of the various 
timescales will lead to temporally induced scatter when- 
ever quantities that depend on these timescales are com- 
pared. 

We plot in Figure [TO] the current central cooling time 
against the average cavity age calculated as the average 
of three standard age estimates: the sound s peed age . 



the buoyan cy age, and refill age (for details, see Rafi^erty 



et al. 2006 1 . For systems with two cavities, an average 
is taken across both cavities. For systems with multiple 
generations of cavities (Perseus and Hydra A), the av- 
erage age is the average difference in ages between the 
inner and outer sets of cavities (i.e. the time between 
outbursts) . We adopt the range in age from the three es- 
timat es as the error. For systems unique to the 'R afferty| 
et al. ([2006) sample, cooling times were derived following 



the procedure described in Section 4.1 



We note that cavities are rarely seen with ages greater 



than ^ 10 yr ( McNamara & Nulsen 2007 1 , perhaps be- 
cause they disrupt on that timescale or fade into the 
background as they rise. Therefore, systems without vis 
ible cavities would likely lie at i nter -outburst times (cav- 



ity ages) above 10 yr in Figure 10 In keeping with this 



conclusion, such systems are shown in Fig[TO]as lower lim- 
its. For three systems in our sample (A496, A1991, and 
NGC 6338), this limit exceeds the central cooling time 
by a significant factor, implying that feedback cannot op- 
erate e ffectively and significant cooling should occur. In 
A1991, |McNamara fc O'Connell] ( [T989| [19921 ) found ev- 
idence for recent star formation (while we do not find a 
positive color gradient in A1991, our images of this sys- 
tem were affected by poor seeing and scattered light, and 
hence the excess blue emission was likely missed). The 
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Fig. 1 0. — The central coo ling time versus the average cavity 
age from |Rafferty et al.|p006| |. The errors in age reflect the range 
in the three age estimates (see text for details). Lower limits are 
shown for systems without visible cavities under the assumption 
that any cavities must have disrupted (or faded below detection 
limits) and are therefore likely older than 10* yr. The line denotes 
equality between the two times. 

two remaining systems do not show evidence of recent 
star formation. It is possible that these systems do have 
cavities, but that they are below our detection limit (in- 
deed, the existing Chandra image of A496 shows some 
evidence of faint cavities) . 

Figure [lO] demonstrates that the average cavity age in 
all systems with detected cavities is less than the current 
cooling time, a result that is consistent with feedback. 
We note that there is evidence of a trend between the two 
times, as expected if systems with shorter cooling times 
need more frequent outbursts to prevent the cooling of 
large amounts of gas. However, the trend is very weak, 
as systems with similar cavity ages have cooling times 
that differ by up to a factor of ~ 100, possibly because 
the true average time between outbursts is not known 
for the majority of systems in our sample (for systems 
without multiple generations of cavities, the current cav- 
ity age is only a lower limit on the average time between 
outbursts). The central cooling time may also vary over 
time as the ICM is heated and cools, inducing additional 
scatter in Figure [TO] 

The timescales relevant for star formation are the time 
for cool gas to form stars and the time for star formation 
to fade below detection limits. The former timescale is 
not well kno wn, but is believed to be on the order of 
10^ 




1.0 10.0 100.0 
Separation (kpc) 



1000.0 



Fig. 11. — Color gradient versus the cooling time at the projected 
radius of the CDG's core. Symbols and the dashed line are the same 
as those in Figure [s] 

objects in our sample) would fade to A{U — I) ^ —0.2 
(typical of the red objects) in a few 10* yr. Alterna- 
tively, the mean change across the sample in [/ — / or 
U ~ R color calculated between the inner and outer radii 
of star formation (instead of between 5-10 kpc or 5-20 
kpc as was done for the A colors) is ~ 0.52, resulting in 
the same fading time. Consequently, detectable emission 
from star formation should typically persist for a few 10* 
yr after cooling ceases. 

The total process of gas condensation, star formation, 
and the fading of the UV emission should therefore re- 
quire ~ 5 X 10*- 10^ yr. This timescale is longer than 
the typical time between outbursts or the typical cen- 
tral cooling time of the systems in our sample, thus 
the colors that we observe could represent the properties 
of star-forming populations averaged over several heat- 
ing events. Additionally, the cooling time at 12 kpc, 
which is comparable to this total timescale, could roughly 
trace the state of the ICM averaged over timescales long 
enough to be typical of the conditions under which the 
stars we see now were accumulated. This scenario might 
explain the strong dependence of star formation on the 
cooling time at 12 kpc seen in Figure [5] 

5.4. Star Formation and the CDG's Location 

It is quite common to find the CDG offset 10 kpc 
or more in projection frort i the X-ray core (e.g., jPatel 



et al. 2006 Edwards et al. 2007). If the CDG does not 



yr (e.g., Mouschovias et al. 2006). The timescale lie at the core of the cooling flow where cooling times 

are short, but instead resides at larger radii where the 
cooling times are relatively long, we should expect little 
active star formation if it is fueled by the cooling ICM. 
To investigate this possibility, we plot in Figure [TT] the 
color gradient against the projected physical separation 
between the CDG's core and the cluster's core. In the 
right panel of Figure l5] we plot the color gradient against 
the cooling time of the ICM at the projected location 
of the CDG's core. In neither plot do the red objects 
segregate to large separations or long cooling times (al- 



for star formation to fade may be roughly estimated 
for our sample using the A colors given in Table [5] and 
the approximately power-law relation between the color 
and age of a young stellar population. For a burst of 
star formation with a metallicity on e-half solar, the stel- 
lar po pulation synthesis models of Bruzual & Chariot 
( 2003 1 predict that the U — I color fades at a rate of 
A[U --I)/Alogt « 1.2 mag dex"-'^. If we assume that the 
average age for the star forming regions in our sample is 
10* yr, a A{U — I) color of 0.3 mag (typical of the blue 
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though those objects with large separations tend to have 
long cooling times). However, the blue objects all have 
small {< 20 kpc) projected separations. We note that 
projection could result in measured separations that are 
much smaller than the true ones. However, on average, 
such effects should be relatively small. 

In summary, although a small (projected) separation 
between the X-ray cusp and the galaxy does not guar- 
antee a blue core, it appears to be a necessa ry condi- 
tion. This result agrees with the findings of Edwards] 



et al. (20071, who find that all the CDGs in their sample 
with strong optical line emission (indicative of an ion- 
izing source, such as young stars or an AGN, near cool 
gas) lie within 50 kpc of the cluster's X-ray core (see also 
Crawford et al."1999|. Therefore, the CDG's location rel- 
ative to the cluster's core is critical to the presence of star 
formation, a finding which supports the hypothesis that 
the cooling ICM, which preferentially cools in the core, 
fuels the star formation. 

5.5. The Suppression of Star Formation by AGN 
Feedback 

The lack of excess blue emission in some objects with 
short cooling times may be due to several factors. In 
general, two situations could apply: either active star 
formation is present, but the excess blue emission is ob- 
scured, or no active star formation is present. In the 
former case, dust obscuration is the most likely cause; 
however, as we argue below, dust is unlikely to obscure 
all signs of significant star formation. In the latter case, 
where no active star formation is occurring, some mech- 
anism for suppressing star formation must be present, 
since cooling times are short and significant condensa- 
tion from the ICM onto the CDG should be occurring. 
An obvious (but perhaps not the only) mechanism for 
the suppression of cooling and star formation in these 
systems is AGN feedback. In this section, we examine 
the red systems with short cooling times and discuss the 
possibility that dust obscuration or suppression by AGN 
feedback could explain their properties. 

Dust is often associated with cold gas and star forma- 
tion and is common in the cores of CDGs (e.g., Laine 
et al. 2003). Dust will preferentially scatter and absorb 



short-wavelength emission, resulting in observed colors 
that are redder than the intrinsic ones. However, dust 
in CDGs is generally obser ved to be patchy or filamen- 
tary (e.g., Laine et al.|2003 ), not spread smoothl y across 
the galaxy in signiticant quantities (see however Silva & 
Wise 19961. Therefore, it is unlikely that dust would 
obscure the entire star forming region if star formation 
in all of the objects with short cooling times is similar. 
New far-infrared measurements of star formation in cool- 
ing flows obtained with the Spitzer o bservatory roughl y 
agree with optical and near-UV rates ( O'Dea et al. 2008 1 . 
Furthermore, for the 17 systems in our samp le that are 
included in recent i nfrared studies ^ Egami et al. 2006 
Quillen et al. 2008!, the 8 systems 



O'Dea et al. 


2008 


that show in 


rared 



rared evidence of star formation also have 
blue cores, implying that star formation in these systems 
is not heavily obscured. Of the remaining 9 systems with- 
out infrared-detected star formation, most (7) have red 
cores. The remaining two, A1991 and MS 1455.0-1-2232, 
have only weakly blue cores. Therefore, near-UV and 
infrared data appear to have similar sensitivity to the 
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Fig. 12. — Color gradient versus the ratio of cavity power to 
cooling luminosity. The vertical, dotted line divides systems with 
cavity power in excess of that needed to balance the cooling lu- 
minosity {Pcav/L-^ > 1) from those with insufficient cavity power 
(fcav/ix ^ 1) to do so. The dashed line and symbols are the same 
as those in Figure [s] 

presence of star formation in these systems, and we are 
unlikely to have missed significant star formation in the 
red systems. 

If instead star formation is suppressed in some systems 
with short cooling times, AGN feedback may be the pri- 
mary cause. As noted in Section |5.2| 42 of the 46 sys- 
tems in our sample have detections at 1400 MHz, and all 
of the systems with short cooling times have detections. 
Therefore, AGNs appear to be almost ubiquitous in the 
cores of cooling flows. Additionally, AGN feedback at 
the cores of clusters has been shown to be energetically 



in many systems ( 


Birzan et al.|2004 Rafferty et al.|2006 


Dunn & Fabian 


2006|). Some systems with recent AGN 



feedback is supplying excess heat above that required 
to balance cooling losses from the ICM, whereas other 
systems appear to be in a coohng stage, in which AGN 
feedback cannot balance the entire cooling luminosity of 
the cluster. Simulations suggest that these systems may 



cycle between the two stages (Omma & Binney 2004 
Ciotti fc: Ostriker|2007 ). An intriguing possibility is that 



the blue systems occupy different temporal locations on 
the cycle of cooling and heating than the red systems. 
The blue systems could reside in clusters in which the 
central AGN is not currently supplying sufficient heat to 
offset cooling, whereas in the red systems, the AGN pre- 
vents large amounts of cooling from occurring, despite 
the sometimes short cooling times. In this scenario, one 
might expect that those objects with large amounts of 
heating relative to cooling would be redder than those 
with insufficient heating. 

To test this prediction, we plot in Figure [T2] the color 
gradient versus the ratio of cavity power to ICM lumi- 
nosity within the coolin g radius, Pcav/^x, w here Pcav 
and Lx were taken from Rafferty et al. ( 2006 ) . Prav was 



calculated as 4py/icav, and Lx is the bolometric lumi- 
nosity of the X-ray-emitting gas within the radius inside 
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which tcooi < 7.7 X 10^ yr. While our sample lacks a large 
number of systems with a high ratio of cavity power to 
cooling luminosity (the X-ray-cavity sample of RafFerty 
|et al. 2006, however, consists of roughly equal numbers 
of systems above and below a ratio of unity) , it appears 
that systems with excess blue emission are more likely to 
have a low ratio of Pcav/Lx (i-e. insufficient heating to 
balance cooling). The average (median) ratio of Pcav/^x 
for objects with detected X-ray cavities and positive gra- 
dients is 0.65 (0.28); for objects with negative gradients, 
it is 3.78 (1.51). However, a K-S test does not rule out 
the possibility that the positive- and negative-gradient 
cavity samples share the same parent distribution of the 
ratio of cavity power to ICM luminosity (the resulting 
probability that they do is 0.17). A larger sample of 
systems with high ratios of Pcav/Lx will be critical to 
test the hypothesis that AGN feedback is quenching star 
formation. 

For comparative purposes, we estimated Pcav for the 
remaining systems (that lack evidence of X-ray cavi- 
ties) with short cooling times (those with icooi[12 kpc] < 
8 X 10* yr) from the monochromatic 1400 MHz r adio lu- 
minos i ties u sing the jet-power scaling relation of IBirzan] 



eFaTlpOMI : 



logPeav = (0.35 ±0.07) log Pi4oo + (1.85 ±0.10), (7) 

where Pcav has units of lO''^ erg and P1400 has units of 
10^'^ W Hz~^. The rad io fluxes for all sys tems were taken 
from the VLA F IRST ( [Becker et al.|l995| or NVSS 
don et al. 1998| catalogs, with the exception of A 2390, 
the hux 01 which was taken from Owen et al. ( 1982 1 . This 
scaling relation was calibrated using a large sample of X- 
ray-cavity systems in clusters (many of which are also in 
our sample) and should therefore be generally applicable 
to our sample. It should be noted that the scatter about 
this relation is large {a w 0.8 dex), and hence the cavity 
power for a given system may be misestimated by a large 
factor. However, averages of the resulting cavity powers 
over many systems should be reasonably reliable. 

Of the 12 systems for which the cavity power was es- 
timated in this way, 10 have evidence of recent star for- 
mation. We include A1991 arn ong these system, since 
McNamara fc 0'ComieIl| ( |1992[ ) found a positive U - I 
gradient in Al99l (6W - Jr ^O.036 ± 0.022), and[Mc' 



Namara & O'Connell 
a weak blue excess; t 



( 1989 1 found spectral evidence of 
lerefore, we have probably missed 
the blue emission in this system (our observations were 
affected by scattered light and poor seeing). These 10 
systems have (Pcav/Lx) = 0.1. The remaining 2 sys- 
tems (A1361 and A496), with negative gradients, have 
(Pcav/Lx) = 0.45. 

Addi tionally, a search of the literature (e.g., Crawford 
fc Fabian.l993i [Roche k Eales|2000l [Hicks fc Mushotzky 



20051 reveals that, o f the 32 cavity systems studied in 
Ratterty et al. ( [2006 1, only 3-6 of the 16 systems with 
Pcav/Lx > 1 show evidence of recent star formation 
whereas 1 0-13 of the 16 syste ms with Pcav/Lx < 1 do 



While the RafFerty et al. ( [2006| sample, which was se- 
lected through visual inspection of publicly available ob- 
servations in the Chandra Data Archive, may be biased 
towards sys tems with easily detected cavities (see e.g., 
Diehl et al. 2008), such biases should not be related to 
the detectability of star formation. Therefore, the ten- 
dency for systems with low ratios of Pcav/Lx to host 



recent star formation and those with high ratios to lack 
such star formation does not appear to be due to the bias 
in our sample towards systems with low ratios. 

The tendency for blue systems to have low ratios of 
Pcav/Lx and red systems to have high ratios lends sup- 
port to the feedback scenario outlined above, and can 
explain 5 of the 7 red systems in Figure [5[ with a cooling 
time at 12 kpc less than - 8 x 10* yr: A1361, A133, 
HCG 62, A2052, MS 0735.6+7421, all of which have 
Pcav/Lx 1- The two remaining systems are A496 
(tcooi[12 kpc] « 7x 10* yr and Pcav/Lx « 0.3) and A2029 
(tcooi[12 kpc] w 6 x 10* yr and Pcav/Lx « 0.07), which 
do not appear to have enough AGN heating to balance 
cooling and yet have short cooling times and are red. It 
is possible that we have underestimated the cavity pow- 
ers for these systems, particularly for A496, for which we 
estimated P^av using equation ([t]). Additionally, 



Clarke 



et al. (20041 found a spiral excess structure in Chandra 
data at the core of A2029 that most likely indicates a 
recent infall or merger, which may have disrupted cool- 
ing and star formation in the core without increasing the 
cooling time at 12 kpc greatly. Lastly, we may have sim- 
ply caught the two systems just before star formation will 
occur (i.e., they are transitioning to a cooling phase). 

We note that the dividing point between blue and red 
systems appears to be roughly at Pcav/Lx ~ 1, implying 
that, if the above scenario is correct, the cavity power 
is a good tracer of the heat input of the AGN. Cavities 
may trace the total heat input well because their powers 
scale with those of associated heating mec hanism s, such 
as sh ocks (e.g 



Nulsen et al. 


2005a [McNamara et al. 


d waves (Sand 


ers & Fabian 


20071 . It 



mechanism in these systems. Nevertheless, it appears 
from Figure [T2| that the assumptions that go into the 
calculation of the cavity power and cooling luminosity 
are approximately correct, or the dividing line would be 
shifted far from unity. The division between blue and 
red systems is not perfectly clean, however, as some blue 
systems fall somewhat above unity and some red systems 
fall below unity. Uncertainties in the cavity timescales 
that go into the calculation of the cavity power and the 
time required for the fading of star formation after cool- 
ing has been quenched would blur the division between 
red and blue systems. 

6. CONCLUSIONS 

We have measured broadband optical colors and de- 
rived X-ray properties for a sample of CDGs in both 
cooling-flow and non-cooling-flow clusters to investigate 
possible connections between the presence of star forma- 
tion and the properties of the cooling flow and the AGN. 
We show that, on similar spatial scales, the presence of 
central blue colors, indicative of active star formation, 
depends critically upon the presence of cooling gas with 
short cooling times. Blue cores are found to occur only 
in the clusters in our sample with central ICM cooling 
times below a threshold of 5 x 10* yr, with central en- 
tropies below ~ 30 kev cm^^, and where the separation 
between the X-ray core and the CDC is less than ^ 20 
kpc. Finally, we find evidence that suggests that the 
lack of observable signatures of star formation in some 
systems with short cooling times could be due to excess 
heating by the central AGN that prevents net cooling 
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from occurring. 

Our results provide compelling evidence for the cycle 
of cooling, star formation, and energetic feedback. This 
cycle may be closely related to the physics responsible for 
several unsolved problems in galaxy formation, including 
the turnover at the bright end of the galaxy luminosity 
function and the preponderance of giant red galaxies, 



rather that giant blue gal axies at late times (Bower et al. 
2006] [Croton et al.|pO%| 
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